Thermal Behavior of Spin Clusters and Interfaces in two-dimensional Ising Model on 

Square Lattice 
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Extensive Monte Carlo study of two-dimensional Ising model is done to investigate the statistical 
behavior of spin clusters and interfaces as a function of temperature, T. We use a tie-breaking 
rule to define interfaces of spin clusters on square lattice with strip geometry and show that such 
definition is consistent with conformal invariant properties of interfaces at critical temperature, Tc. 
The effective fractal dimensions of spin clusters and interfaces (dc and dj, respectively) are obtained 
as a function of temperature. We find that the effective fractal dimension of the spin clusters behaves 
almost linearly with temperature in three different regimes. It is also found that the effective fractal 
dimension of the interfaces undergoes a sharp crossover around Tc, between values 1 and 1.75 at 
low and high temperatures, respectively. We also check the finite-size scaling hypothesis for the 
percolation probability and the average mass of the largest spin-cluster in a good agreement with 
the theoretical predictions. 
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Two dimensional (2D) Ising model as a solvable pre- 
scription model in hand, and its extension to g-state 
Potts model [1] have been the subject of intense research 
interest for decades. Many of their thcrmodynamical pa- 
rameters and behaviors can be characterized in terms 
of some fractal geometrical objects, e.g., spin clusters 
and domain walls. Most of studies have been focused 
to describe the behavior of these models at critical tem- 
perature Tc, at which they exhibit a continuous phase 
transition (for (? < 4), and less attention is made to in- 
vestigate off-critical characterization of such systems at 
temperatures far from Tc. At T = Tc, conformal field 
theory (CFT) plays an important robust role to describe 
the universal critical properties in two dimensions. Be- 
sides CFT, theory of stochastic Loewner evolution (SLE) 
invented by Schramm [2] provides a geometrical under- 
standing of criticality which states that the statistics of 
well defined domain walls (or curves, e.g., spin cluster 
boundaries in 2D Ising model) in upper half plane EI is 
governed by onc-dimensional Brownian motion (to re- 
view SLE, see [3]). Therefore it is expected that for 
example in 2D Ising model, the geometrical exponents 
such as the fractal dimension of a spin cluster and its 
boundary as well would be related to the thcrmodynam- 
ical exponents [4] . The study of the fractal structure and 
the scaling properties of the various geometrical features 
of the Ising model has been subject of huge scientific 
literature (see for example [5, 6, 7, 8, 9] and references 
therein) . It is also well known that most two-dimensional 
critical models renormalize onto a Gaussian free field the- 
ory (Coulomb gas) [10]. Many exact critical exponents 
have been computed by using the Coulomb gas technique 
[11]. These include various geometrical exponents of two- 
dimensional Ising model [12], and general g-state Potts 



model [13, 14]. 

The geometrical objects reflect directly the status of the 
system in question under changing the controller param- 
eters. Temperature can play the role of such controller 
parameter in 2D Ising model. 

Investigation of the dependence of geometrical expo- 
nents in 2D Ising model, equivalent to q = 2 states Potts 
model, is the main subject of the present paper. To be 
consistent with the postulates of SLE at T = Tc, we con- 
sider the model on a strip of size x Ly, where is 
taken to be much larger than Ly = L, i.e., = 8L. We 
simulate the spin configurations of 2D Ising model on 
square lattice using Wolf's Monte Carlo algorithm [15], 
based on single cluster update. Before going into the fur- 
ther details, let us address an ambiguity that arises when 
one intends to define an Ising interface on a square lat- 
tice, and then introduce a rule which seems to produce 
well-defined interfaces on square lattice. The importance 
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FIG. 1: An Ising interface defined on hexagonal lattice cor- 
responding to a spin configuration on triangular lattice with 
a fixed boundary condition at the real line in H, as explained 
in the text. 
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of such definition backs to its relevance to the SLE in- 
terfaces at criticahty. As will be discussed later, it is 
believed that the critical Ising interfaces can be defined 
by the theory of SLE in the scaling limit [16, 17] . Thus we 
need to have a unique procedure to define operationally 
the hulls of the Ising spin clusters without any self inter- 
section and ambiguity. However we will use our following 
procedure to define well defined interfaces in Ising spin 
model, one can simply extend it for any two-dimensional 
model defined on square lattice, e.g., for interfaces of 
general g-state Potts model or contour lines of random 
growth surfaces [18] etc, with appropriate substitutions 
of spins up and down. 

Consider an Ising model on a triangular lattice in up- 
per half plane on which each spin lies at the center of 
a hexagon having six nearest neighbors and the spin 
boundaries (defining the interface) lie on the edges of the 
honeycomb lattice (see Fig. 1). To impose an interface 
(which separates the spins of opposite magnetization), 
growing from the origin on the real line to infinity, a 
fixed boundary condition can be considered in which all 
spins in the right and left sides of the origin are up ('+') 
and down ('—'), respectively. The Gibbs distribution in- 
duces a measure on these interfaces. 
To define an interface, a walker moves on the edges of 
the hexagonal lattice starting from origin at the bottom. 
At each step the walker moves according to the follow- 
ing rule: turns left or right according to the value of the 
spin in front of it ('-|-' or '— ', respectively). The resulting 
interface is a unique interface which never crosses itself 
and never gets trapped. Such a interface, at T = Tc, is 
believed to be described by SLE in the continuum limit 
[16, 17]. 

This procedure to define the interface should be mod- 
ified for spin configuration on square lattice. This is be- 
cause that there are some choices for the square lattice, 
at places with four alternating spins. We first introduce 
a tie-breaking rule which the walker regards at each step 
and then we show that this definition is consistent with 
the predictions of SLE for such interfaces at T = Tc- 

Consider a spin configuration on a strip of square lat- 
tice in H, with the same boundary conditions as above. 
A walker moves along the edges of the dual lattice (the 
lattice shown by the dotted-dashed fines in the Fig. 2), 
starting from the origin. According to the boundary con- 
ditions at the first step of the walk that the spins lie 
at the right of the walker, this direction is chosen to be 
the preferable direction. After arriving to each site on the 
dual lattice, there are three possibilities for the walker: 
it can cross one of the three nearest bonds of the original 
lattice. At the first step of selection, it chooses the bonds 
containing two different spins that crossing each of which 
leaves the spin at the right and spin '— ' at the left of 
the walker. The directions right and left are defined lo- 
cally according to the orientation of the walker. After the 
first selection, if there are yet two possibilities to cross, 
the walker chooses the bond which accords with the turn- 
right tie breaking rule: it turns towards the bond which is 
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FIG. 2: An Ising interface defined on square lattice, dual of 
the original square lattice including a spin configuration, with 
a fixed boundary condition at the real line in H. The interface 
is generated applying the turn-right tie-breaking rule. The 
same procedure can be used to define such interface for down 
spins (' — ') according to turn-left tie-breaking rule. 



in its right hand side with respect to its last direction at 
the last walk; if there is not any selected bond at its right, 
it prefers to move straightly and if there is not also any, 
it turns to its left. The procedure is repeated iteratively 
until the walker touches the upper boundary. The result- 
ing interface is again an interface which touches itself yet 
never crosses itself and never gets trapped. The same 
procedure can be used to define another interface with 
Ze/t-preferable direction as turn-left tie-breaking rule. 
It would be worth to mention that the procedure in- 
troduced here yields not just a unique cluster boundary 
without any ambiguity on the square lattice, but one can 
check that, any other definition for the interface leads to 
an incorrect boundary of the cluster (for example at ver- 
tices with more than one possibility, just these introduced 
options lead to the 'true' boundary of the considered clus- 
ter and any other option, for example choosing randomly 
the directions left or right, may enter the boundary of a 
spin which does not belong to the cluster. Note that a 
spin cluster is defined as a set of nearest neighbor con- 
nected sites of like sign.). 

Let us now show that the resulting interface is compati- 
ble with the properties which comes from their conformal 
invariant nature at T = Tc- 

The Wolf's Monte Carlo algorithm is used to simulate 
the spin configurations at T = Tc, on the strip of square 
lattice and of aspect ratio 8, and boundary conditions as 
discussed above. For each size L, about Monte Carlo 
sweeps are used for equilibration. An ensemble of 2 x 10^ 
independent samples is collected for each sample size L, 
where each of which was taken after lOL Monte Carlo 
steps. 

Each spin cluster has been identified as a set of con- 
nected sites of the same spin using Hoshen-Kopelman 
algorithm. We just take the samples including a vertical 
spanning cluster in the y-direction. Then an ensemble 
of corresponding spanning interfaces was obtained using 
mentioned turn-right (left) tie-breaking rule. 
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The fractal dimension of the interfaces at this critical 
temperature, dj(Tc), is obtained using the standard finite 
size scaling. The length of an interface I scales with the 
sample size as / ~ L'^dTa) ^ ^j^g fractal dimension of con- 
formally invariant curves is provided by SLE [3] generally 
as d/ = I + k/8, where diffusivity k classifies different uni- 
versality classes, and for Ising spin-cluster boundaries it 
is conjectured to be k = 3 and thus dj{Tc) = — 1.375. 
As shown in Fig. 3, the best fit to our data collected 
for sizes 30 < L < 500 yields the fractal dimension 
d/(Tc) = 1.371 ± 0.005. 

Another prediction of the theory of SLE for such criti- 
cal interfaces is the winding angle statistics [2] . We define 
the winding angle 9 as defined by Wieland and Wilson 
[19]. For each interface we attribute an arbitrary wind- 
ing angle to the first edge (that wc take zero). Then the 
winding angle for the next edge is defined as the sum of 
the winding angle of the present edge and the turning 
angle to the new edge measured in radians. It is shown 
that [19, 20] the variance in the winding grows with the 
sample size like 

(0=a+JlnL, (1) 

where k = %[dj[Tc) — 1], and a is a constant whose value 
is irrelevant. So the exact value of k for critical interfaces 
of 2D Ising model should be k = 3. 

The figure 4 indicates that our result for k is in a good 
agreement with the predicted value. Wc find that n = 
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FIG. 3: (Color online) Log-log plot of the average length of 
a spanning interface generated using the tie-breaking rule 
introduced in the text, versus the wide of the strip L, at 
critical points. Main: for Ising model. Inset: for the hull 
(the upper graph) and its external perimeter (EP— the lower 
graph) of critical site percolation. The values of the best fit 
to the data are represented aside each one, with an error of 
~ 0.005. 
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FIG. 4: (Color online) Variance of the winding angle for span- 
ning interfaces generated using the tie-breaking rule intro- 
duced in the text. The solid line is set according to the Eq. 
(1), with a = —0.29 and k — 3. In the inset, the variance in 
semilogarithmic coordinates. 



3.012 ± 0.005. 

We have also tested other conformal invariant properties 
of the interfaces such as Schramm's formula for the left 
passage probability of the interfaces, consistent with the 
theory (the results are not shown here). 

To investigate another concern about the systems with 
more complicated interfaces, we did such experiments for 
the critical site percolation [21]. The fractal dimension 
of the hull and its external perimeter are obtained as 
df = 1.751±0.002,and df^ = 1.335±0.002, respectively 
(see Fig. 3) in a good agreement with the duality relation 
predicted from the conformal invariant property [22] 

(df-i)(dr-i) = i. (2) 

In the rest of the paper, let us consider the statistical 
geometrical response of the Ising model to the tempera- 
ture. We show experimentally that how the statistics of 
the spin clusters and their boundaries behave as a func- 
tion of temperature. We try to measure the correspond- 
ing fractal dimensions at length scales smaller than the 
correlation length ^, using the standard finite size scaling 
as done at critical temperature above. 
Features of the spin clusters at three different tempera- 
tures are shown in Fig. 5. These represent what we ex- 
pect to happen: at zero temperature, because of the used 
boundary conditions, the ground state of the spin con- 
figuration splits the system into two segments, one with 
spins up and the other with spins down which are sepa- 
rated with a straight interface. Increasing in the temper- 
ature induces a fractal random feature on spin clusters 
and interfaces. The interfaces are some non-intersecting 
curves (in H) which can be described, in the continuum 
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FIG. 5: (Color online) Spin clusters of 2D Ising model on strip of square lattice with size of L = 120, and aspect ratio 6, at 
different temperatures from top to bottom: T — Tc = —0.2, and, 0.2. The boundary conditions (be) used for simulation are 
fixed for the lower boundary, antiperiodic at sides and free be for upper one. The spin-down clusters are shown white. The 
be imposes an interface at the boundary of the spanning cluster (dark colored) starting from the origin (using the turn-left 
tie-breaking rule in these figures) and ending at the upper boundary. As temperature increases the interface gets more space 
filling. 
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FIG. 6: (Color online) Main: Log-log plot of the average 
length of a spanning interface I, versus the wide of the strip L, 
at two different temperatures T — Tc = —0.1 and -1-0.1. These 
graphs show the scaling property and the fractal behavior 
of the interfaces far from criticality at length scales less than 
the correlation length. Inset: Log-log plot of the average mass 
of a spanning cluster M, versus the wide of the strip L, at 
T - Tc = -0.1 and -fO.l. The graph for T - Tc ^ +0.1 is 
shifted upwards by 2. 



limit, by a dynamical process called Loewner evolution 



[23] with a suitable continuous driving function Q as 
dgtjz) ^ 2 

where, if we consider the hull Kt, the union of the curve 
and the set of points which can not be reached from in- 
finity without intersecting the curve, then gt{z) is an an- 
alytic function which maps M \ Kt into the EI itself. 
At zero temperature the driving function (^t , is an specific 
constant, at T = Tc it should be proportional to a stan- 
dard Brownian motion Bt as Ct = ^/kBi with k = 3, and 
it may be complicated random function at other different 
temperatures. 

At high-tcmpcrature limit, each spin gets the directions 
up or down with probability p = \/2 and so, it is con- 
jectured to correspond to the critical site percolation on 
triangular lattice (on which the percolation threshold is 
exactly at pc = 1/2), and it is expected that the driving 
function converges to a Brownian motion with diffusiv- 
ity of K = 6. For the case of square lattice, since the 
percolation threshold in two dimensions is at Pc ^ 0.59, 
so at high-tcmpcrature limit where p = ^ < Pc, the sys- 
tem will be below the threshold and the crossover to the 
critical site percolation will not be seen any more. 

Before looking at the temperature dependence of the 
fractal dimension of the spin-clusters, let us discuss more 
about their scaling properties from the point of view of 
theoretical expectations. 

The Ising model is expected to be scale-invariant (on 
scales much larger than lattice spacing a) only at renor- 
malisation group fixed points, i.e., T = Tc and T = oo 
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(on triangular lattice) . At those points one expects well- 
defined power-law behavior for clusters and their hulls 
on all scales L ^ a. For T just above Tc, where the cor- 
relation length ^ is finite and ^ 3> a, one expects to see 
behavior characteristic of the critical point Tc on scales 
a ^ L <C ^, and of the high-temperature fixed point on 
scales L ^ 

Thus, according to the theory, there should be no such 
thing as Hhe fractal dimension at temperature T\ except 
for T — Tc and T = oo , instead one should see a crossover 
between two different values. If one chooses a sufficiently 
narrow range of length scales one will see an effective 
fractal dimension, which will have the appearance of de- 
pending on temperature. However, for the Ising model on 
square lattice, since the crossover to the critical percola- 
tion at high-tcmpcratures no longer exists, the behavior 
of the effective fractal dimensions is governed by just the 
behavior at T = Tc for length scales a ^ i ^ ^. In 
order to determine the behavior of such effective frac- 
tal dimensions as a function of temperature, we measure 
them in an almost narrow range of sizes L, which seem to 
be much smaller than the correlation length and within 
the range the scaling properties are held. 

Figure 6 shows the procedure we perform to measure 
the effective fractal dimension of the spin clusters and in- 
terfaces at different temperatures. The finite size scaling 
reduces substantially the statistical errors in estimating 
the fractal dimensions. The average is taken over 10''^ in- 
dependent samples of aspect ratio 4. at each temperature 
below Tc for each sample size (only the spanning clus- 
ter in each configuration and the corresponding interface 
was considered). Since the probability to have a span- 
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FIG. 7: (Color online) Effective fractal dimension of spin clus- 
ters as a function of temperature. It changes almost linearly 
in three difterent regimes: low temperature with dimension of 
2, rapid decreasing around Tc and a crossover to a different 
linear behavior far from Tc. The slope of the dashed- lines 
differs by one order of magnitude. Each point is obtained 
using finite size scaling for 10 different sizes in the range of 
50 < L < 500. The error is less than the symbol size. 
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FIG. 8: (Color online) Effective fractal dimension of spin clus- 
ter boundaries as a function of temperature. The error is less 
than the symbol size. 



ning cluster diminishes when temperature increases (as 
will be discussed later), the average is taken over 2 x lO"*^ 
independent samples for T > T^ and the samples were 
gathered on strip of aspect ratio 8. 

The exact values for the fractal dimensions of spin clus- 
ters and interfaces are known just for at critical tem- 
perature Tc, as dc{Tc) = ^ = 1.9479... and di{Tc) = 
^ = 1.375, respectively. Our measurements of fractal 
dimensions at Tc which give dc{Tc) = 1.9469 ± 0.001 
and di{Tc) — 1.371 ± 0.005 are in a good agreement 
with the exact results. These values were obtained for 
30 < L < 500. The same measurements for T ^ Tc, were 
done for 10 different sizes within 50 < L < 500 (the ex- 
amples are shown in Fig. 6). 

Fig. 6 shows the scaling properties and the fractal behav- 
ior of the spin clusters and interfaces at T ^ Tc, within 
the selected range of size. 

To quantify the geometrical changes of the spin clusters 
at different temperatures, we measure the effective fractal 
dimensions of the spin clusters and their perimeters. At 
each temperature, we use the scaling relation between 
the average mass of the spanning spin-cluster M , and the 
width of the strip L, to measure the fractal dimension of 
the spin-clusters — i.e., M ~ L'^". 

Corresponding fractal dimension of spin clusters as a 
function of temperature is shown in Fig. 7. This suggests 
three different regimes, one for low temperatures in which 
the dimension of the spin clusters is 2. The second regime 
is in the vicinity of the critical temperature: a linear 
dependence of the fractal dimension on temperature with 
a sharp decreasing which is governed by criticality. A 
crossover happens at temperature above critical region 
which changes the slope of the linear decrease by about 
one order of magnitude at high temperatures. 

Such a crossover can be also seen in the behavior of the 
effective fractal dimension of the interfaces as a function 
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of temperature. As shown in Fig. 8, at low temperatures 
the effective fractal dimension of the interfaces is close 
to 1 and it increases with temperature. In the vicinity 
of the critical temperature it increases again sharply and 
then crosses over to the value very close to 1.75, which 
is the fractal dimension of the hull of critical percola- 
tion. The whole behavior looks like a hyperbolic tangent 
function. The other theoretical predictions for the ge- 
ometrical features considered in this paper and we are 
interested in checking them, are about the percolation 
observables. The finite-size scaling hypothesis states that 
the percolation probability Pg i.e., the probability to have 
a spanning cluster at temperature T, reaching from one 
boundary to the opposite one, behaves like [6] 

= P.(£/0, (4) 

where the correlation length behaves like ^ ~ (T — Tc)"", 
with ly = 15/8 for the Ising spin geometric clusters. 
In order to investigate this hypothesis, we have done sim- 
ulations of Ising model on square lattices of different size 
with free boundary condition, and the measurements 
are taken by averaging over 2 x 10^ independent samples 
at each temperature. As shown in Fig. 9, curves mea- 
sured on lattices of different size all cross at the critical 
point (in the figure this observable is shown as a function 
of the inverse temperature (3). As can be seen from the 
figure, applying the scaling theory Eq. 4, results data 
collapse onto a single function, in a good agreement with 
the theoretical predictions. 

The other observable we consider is the scaling behav- 
ior of the average mass of the largest spin-cluster, M. 
According to theory, this should have the scaling form 

M ^ L'^-(^=)p(P/^), (5) 




FIG. 9: (Color online) Finite-size scaling plots of the data 
for the percolation probability, measured on square lattices 
of different size . Inset: the percolation probability as a 
function of the inverse temperature /3. 




FIG. 10: (Color online) Data collapse for the average mass 
of the largest spin-cluster Af, measured on square lattices of 
different size . Inset: the strength of the largest spin-cluster 
as a function of the inverse temperature /3. 



where the scaling function F{x) goes to a constant as 
X -> (at T = Tc). 

The suitably rescaled mass of the largest spin-cluster as 
a function of the reduced inverse temperature is plotted 
in Fig. 10, implying the data collapse onto a universal 
curve. 

In conclusion, we studied the geometrical changes of 
the spin clusters and interfaces of two-dimensional Ising 
model on square lattice in the absence of external mag- 
netic field, as a function of temperature. We introduced a 
well-defined tie-breaking rule to generate nonintersecting 
interfaces on square lattice, which are shown to be consis- 
tent with the predictions of conformal invariance at the 
critical point. The results are also checked for critical 
site percolation in a good agreement with the analytical 
predictions. 

We also investigated the effect of the temperature on the 
statistical properties of geometrical objects by measuring 
the effective fractal dimensions of the spin clusters and 
interfaces as a function of temperature. We showed that 
a crossover happens which distinguishes between the be- 
havior of these geometrical objects near the critical tem- 
perature and that of at high temperatures. 
We also applied the finite-size scaling hypothesis for both 
the percolation probability and the average mass of the 
largest spin-cluster, and we found a data collapse onto a 
universal curve, in a good agreement with the theoretical 
predictions. 

I am indebted to J. Cardy and D.B. Wilson for elec- 
tronic discussions and their constructive comments. I 
also would like to thank H. Dashti-Naserabadi for his 
helps on programming. 
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